﻿ Esantionarea semnalelor http://shannon etc upt ro/teaching/ps/Cap9 1 pdf http://shannon etc upt ro/teaching/ps/Cap9 2 pdf •Discretizarea variatiei in timp a semnalului, numita esantionare •Semnale de banda limitata •Problema reconstruirii semnalulelor analogice din semnalul esantionat •Teorema esantionarii: esantionarea ideala 1 Teorema esantionarii 1        u t t t              22      Un esantion al lui x(t) este obtinut prin produsul semnalului analogic x(t) cu un impuls foarte scurt dreptunghiular u(t ) de arie 1: Δ 0x t u t x u t          Un alt esantion poate fi obtinut daca se plaseaza impulsul la momentul de timp kT s      s s sx t u t kT x kT u t kT 2 1 Procedura de esantionare: esantionam semnalul cu un tren de impulsuri foarte scurte              s s sx t u t kT x kT u t kT kk      Pentru valori mici ale lui Δ: Distributia Dirac limu t t      0  limu t kT t kT x kT t                0sTs s s kk      Distributia Dirac periodica Esantionarea ideala a semnalului x(t) :  ˆ x t x t t x kT t kT               sTss k   3 Esantionarea ideala Modelul matematic  x t x t t x kT t kT               sT s s k   Modelul esantionarii ideale 4 2 Spectrul semnalului esantionat ideal Semnalul esantionat ideal:  x t x t t x kT t kT               sT s s k   Spectrul semnalului esantionat ideal: ˆ X x t t           sT 5 Aplicand teorema produsului :  1 2 2  ˆ X x t t X k                    sT 2TT kss        1 2 1 2  ˆ X X k X k                   kks s s sT T T T          Spectrul semnalului esantionat apare prin insumarea unor lobi spectrali indexati dupa k 6 3           ssx t x kT t kT k    12  ˆ X X k        kssTT    Spectrul semnalului esantionat ideal este repetitia periodica a spectrului semnalului original Perioada este invers proportionala cu pasul de esantionare T s 7 Din spectrul semnalului esantionat ideal nu se mai poate recupera spectrul semnalului original, deoarece apare o “amestecare” spectrala (aliere) Semnalul original trebuie sa fie de banda limitata pentru ca aceste erori sa nu apara Spectrul semnalului original Spectrul distributiei Dirac periodice Spectrul semnalului esantionat ideal    12ˆ8 X X k       Eroare de aliere kssTT   4 Teorema esantionarii semnalelor de banda limitata x(t)-de banda limitata X0 daca      M  9 •Lobii spectrali successivi nu se suprapun si spectrul semnalului original poate fi recuperat prin FTJ din spectrul semnalului esantionat daca  2  sM 10 5 FTJ ideal tsin c           thpH c t Utilizam un filtru trece jos cu frecventa de taiere ω si amplificarea in banda T cs Eroarea de aliere poate fi evitata Frecventa de esantionare 2  sM  Frecventa de taiere pentru FTJ        csMM 0HT  Conditia de reconstructie perfecta:  rs 11 Raspunsul in frecventa al filtrului de reconstructie este: T,    cs                 csMM rsH T p c0,   c   Semnalul reconstruit este : ˆ       rrx t x t h t cu spectrul: ˆ           rrX X H  1 X k T p X              css ksT   Spectrele egale implica egalitatea aproape peste tot a semnalului reconstruit cu semnalul initial , x t x t a e w  12 r 6   Daca nu se respecta conditia de esantionare , apar2Ms erorile de aliere      MMs13 Reconstructie sint c H T p h t T           cr s r s t ˆ  rrx t h t x t sint c      s s sT x kT t kT  kt sint cx kT T t kT     s s s  kt sint kT  cs   ssx kT T   kst kT sin2t kT  ccs x kT  s     kscst kT 14 7 Frecventa de esantionare minima este  2 si poarta  sM  denumirea de frecventa de esantionare Nyquist In cazul esantionarii la frecventa Nyquist formula de reconstructie devine: sint kT   M s      rsx t x kT   k   Mst kT 15 Teorema WKS (Whittaker, Kotelnikov, Shannon) Daca semnalul este de banda limitata la ,in sensul ca 0x t X    M pentru , atunci xt   este unic determinat de multimea esantioanelor  M   frecventa de esantionare este cel   sMx nT n Z , daca ssale 2 , adica dublul frecventei maxime In conditiile de mai sus semnalul initial xtputin  se poate reconstitui din esantioanele sale, a p t prin relatia: sin2t kT  csc   sx t x kT     ks c st kT cu conditia ca astfel sa fie ales incat sa satisfaca relatia :        cM c s M   16 Nyquist: sFrecventa 2M 8 Reconstructia prin filtrare trece- jos ideala Semnalul se poate reconstrui din curbe de tipul sin x / x •Intr-un punct de esantionare, kTs, suma se reduce la x(kTs) •Intre punctele de esantionare, reconstructia se obtine prin aportul tuturor termenilor sumei 17 •Operatia de reconstructie se mai numeste si interpolare •In punctele de esantionare, un singur esantion defineste semnalul sinT n k  Ms       ssx nT x kT    kMsT n k   sinnk    s                ,ssMssnkx nT x kT x kT x nT 2kknk         1, pentru n k   nk , 0, pentru n k  •Esantionare cu frecventa Nyquist •Intre punctele de esantionare, reconstructia se 18 obtine prin aportul tuturor termenilor sumei 9 Reconstructie prin interpolare liniara •Este posibila si o reconstructie aproximativa a semnalului, prin unirea punctelor determinate de valorile esantioanelor cu linii drepte Semnalul reconstruit x(t) este doar aproximativ r egal cu x(t) •Raspunsul la impuls al filtrului de reconstructie: triunghiular 19 2 2    sinsT sin  2s      r s sH T T  sT   2 s  ideal Raspunsul in frecventa al filtrului de reconstructie 20 10 Reconstructia prin extrapolare de ordinul zero •Este posibila reconstruirea prin extrapolare de ordin zero: valoarea esantionului curent se mentione pana la aparitia noului esantion TT 2sin sinTTssss    2222jjsT       Trsh t p t e e T  s    22sT 2   sinsT sinTs j   22jss H e T e     rsT21 s  2 s •Diferenta dintre filtre (cel de reconstructie si ideal) este semnificativa  erori de reconstructie mari Modulul raspunsului in frecventa al filtrului de reconstructie 22 11 Spectrul semnalului reconstruit:  sin   j  ss                 rrsX X H e X k  k    s 23 Lobii spectrali sunt puternic deformati Chiar daca dupa reconstructia de ordin zero am aplica un filtru suplimentar, trece jos ideal, tot nu s-ar putea recupera semnalul initial Evident pentru ω>>2 ω, erorile din lobul central pot fi mult sM diminuate 24 12 Esantionarea ideala a semnalelor periodice •Consideram semnale periodice de banda limitata la a N-a armonica Cea mai mare frecventa din spectrul lor este 2   N;   00M T 0 •Se esantioneaza semnalul cu o frecventa   MM, s  0 •Prin esantionare lobii spectrali se repeta •Pentru ca sa nu apara suprapunerea lobilor spectrali : N N M N          0 0 0s 25 N   3  00        0 0 0min2 3 7s NN       00cs 26 13 N N M N          0 0 0s     dintre 00Diferenta si trebuie sa fie de forma:M N N         0 0 0, 1,2,M N N R R             0 0 02 2 , 1, 2,sMM N R R R Frecventa minima de esantionare  2 1 2N         00minMs27 Reconstructie prin filtrare trece jos ideala T,    sc     ;H T p      rs c0,   c   NN       00cs Pentru a evita aparitia erorilor de aliere este necesar ca:     00sNN 22N    sM  0   spre deosebire de semnalele aperiodice unde 2s M Pe perioada celei mai rapide componente spectrale, trebuie sa prelevam mai mult de doua esantioane din spectru Daca esantionarea mentine periodicitatea atunci pe perioada celei mai rapide componente din spectru trebuie sa prelevam 3 esantioane 28 (cel putin) 14 T - perioada fundamentalei 0 Esantionarea se face cu ω=(2N+R)ω atunci: s0 22 2NR   0sTT 0T ; 1, 2,TR s 2NR Doar 2N+R esantioane pot fi distincte ca urmare a periodicitatii semnalului supus esantionarii Toate pot fi prelevate intr-o singura perioada a fundamentalei T 0 29 Principiul osciloscopului Acelasi rezultat se poate obtine prelevand esantioanele succesive din perioade diferite          00s s sx kT x T kT x KT kT 1     0 0 0ssT KT kT KT T 2NR Aceasta posibilitate este valorificata in constructia osciloscoapelor cu esantionare 30 15  2N  http://www jhu edu/~signals/sampling 0s NN       00cs 100 /rad s s 50 /rad s c 120 /rad s s 31 Relatii energetice Pentru semnale aperiodice esantionate, este valabila relatia de tip Rayleigh  22      ssW x t dt T x kT  k    Pentru semnale periodice esantionate, relatia de tip Parseval 1M 2211 ; 2 ; 1, 2,P x t dt x kT M N R R         s   0TTM 00k Energia sau puterea pot fi calculate fie din forma de variatie in timp, fie in frecventa 32 16 Esantionarea cu memorare ˆ x t x t t h t x t h t                  sT  tt   2sin sintt    jjt  2222          th t p t e e t  t  22 2 a) Semnal esantionat b) Model matematic al sistemului de modulare a impulsurilor in amplitudine (PAM) 33 Valoarea esantionului prelevat se pastreaza pentru un interval de timp Δt≤Ts Spectrul semnalului esantionat cu memorare ˆ x t x t t h t x t h t                  sT Spectrul semnalului esantionat cu memorare se poate scrie:  tt   sin sintt     jjt 22122                 ssX e t X k e X k     kkssttTT      22 t sint sint jTtt 22s            ssX e X k X k   kkssttTT   2T s 34 17 Lobii spectrului semnalului esantionat cu memorare sunt deformati, 2  nu se poate reconstrui semnalul initial prin FTJ ideala Pentru ca  M lobul central sa fie putin afectat este necesar ca: t  Scurtarea duratei impulsurilor, Δt Reconstructia prin extrapolare de ordin zero este un caz particular PAM, cu 35 Δt=Ts Esantionarea naturala Esantionare cu “decuparea” unor portiuni din semnal; se aplica in multiplexarea in timp a semnalelor analogice time division multiplexing (TDM) 36 18 x t x t q t x t h t t                 ssTT  x t h t kT x t h t kT            ss  kk      •Semnalul este inmultit cu un tren de impulsuri periodice q(t) Ts t 2sint jt  22             th t p t H e  22 37 Spectrul semnalului esantionat natural X x t h t t      sT t  sint j  2122 X e t k           s  2ktT s  2 kt sinskt s jt 22 X e k         s  skktT s 2 factor numeric kt sinskt s jt 22 eXk     s sktT  sk 38 2 factor numeric 19 sT    c, Ht     •Lobul central este asemanator cu cel obtinut 0,   c  prin esantionarea ideala, (nu apar distorsiuni ca in cazul esantionarii cu memorare)       •Reconstruire prin FTJ ideala M c s M 39 Relatia dintre spectrul unui semnal discret si spectrul semnalului analogic din care provine Semnal analogic si spectrul sau Semnal discret si spectrul sau   x t X x n x nT X;                a a a sdd Doua expresii echivalente pentru spectrul semnalului analogic:  12 ˆ ˆ          x t x t t X X k,             sa s sT kssTT     ˆ ˆ x t x t t X x t t kT          sasT k    x t t kT x kT t kT         a s a s s kk    x kT t kT     a s s k  sj kT40 x kT e  as k 20  1j s kT ˆ             a s a sX X k x kT e kksT      Spectrul semnalului discret este:  j k j k             asddX x k e x kT e kk      Se observa ca:     sj kT jk       a s a sx kT e x kT e  kk      sT Relatia dintre spectrul unui semnal discret si spectrul semnalului analogic din care provine este :  12        adX X k    ks s sT T T41 Intre cele doua axe de frecventa corespunzatoare spectrului esantionat respectiv semnalului discret exista relatia:   sT  Se aplica si acum natura periodica a spectrului semnalului discret  dX Intre frecventele maxime din spectru exista relatia:  TT   ; ssMM 42 M 21 Esantionarea semnalelor discrete Problema: In prelucrarea numerica a semnalelor apar situatii in care, ulterior achizitionarii esantioanelor, se constata ca frecventa de esantionare a fost prea mare Solutie: In astfel de situatii, cand nu se mai poate esantiona semnalul analogic, este posibila esantionarea semnalului numeric, retinandu-se tot a N-a valoare Se reduce astfel numarul de esantioane dupa ce semnalul a fost deja esantionat = esantionarea semnalului discret Semnalul discret esantionat se obtine prin produsul:  ˆ                Nx n x n n x n n kN k    x kN n kN        43 k   Pas de esantionare N=3 44 22 Spectrul semnalului esantionat  2             nk;          sNs s s s kN   ˆ           NX x n n  1 2 2 Xk ;             ss  2NN k    Spectrul semnalului discret esantionat este periodic, de perioada Ω s 11NN 1 1 2 ˆ           X X k X k,            r r rs s s  00kkN N N Restrictia pe 2π (perioada spectrului semnalului x[n]) 45 Pas de esantionare N=3 X Spectrul  periodic, de perioada 2π  ˆ X Spectrul  periodic, de perioada Ω s 46 23 Semnalul analogic, x(t), de banda limitata ω , a fost esantionat ideal, cu M pasul de esantionare T  semnal numeric de banda limitata x[n], care s este si el esantionat cu pasul N  semnal numeric esantionat  ˆ xn  Se dovedeste ca semnalul x[n] este esantionabil cu pasul N adica lobii spectrali nu se suprapun, daca semnalul analogic x(t) ar fi putut fi esantionat si cu perioada T’=NT respectandu-se teorema WKS Initial ss semnalul analogic a fost supra-esantionat Frecventa maxima este Ω=ωT, Ω=2π/N, atunci conditia de MMsS reconstruire perfecta este: Ω≥ 2Ω sM  47   s s s sNT T T NT; MM Eroarea de aliere      MMs 48 24 Reconstruirea semnalului discret din esantioanele sale Daca lobii nu se suprapun, prin filtrare trece-jos ideala:  2N , k    c H           r M c s M  0 in, rest 49  Caracteristica in frecventa a FTJ ideal in timp discret (capitolul Filtrare) 1, 2 , sinkkn       cc  H h n           0, in rest n  50 25 Esantionarea si decimarea unui semnal discret ExDupa esantionarea unui semnal discret, intre emplu, N=2 doua valori retinute (esantioane) sunt intercalate N-1 zerouri, care nu aduc nici o informatie despre semnalul esantionat Ele pot fi omise Esantionarea semnalului discret ˆ  Rezulta semnalul decimat  Dxn Relatia dintre semnalul discret decimat si cel esantionat : ˆ     Nnxnx ˆ D Decimare Din semnalul decimat se poate reconstrui semnalul nedecimat prin inserarea N-1 zerouri intre doua valori consecutive 51 Spectrul semnalul decimat :     jm      Nj n j n ˆˆ ˆ ˆ ˆˆ            DDDx n x n e x nXXN e x m e     n n mN 1N 1 Spectrul semnalul digital Spectrul semnalul digital ˆ          sX X k esantionat initial X(Ω) 0kN 1N 12k    ˆ    rDrXX   0kNN  N  0 : 1/ / 0Mk N X N           Intinderea lobilor spectrali ai semnalului decimat este de N ori mai mare decat intinderea lobilor spectrali ai semnalului nedecimat x[n] 52 26 exemplu N=2 Spectrul semnalului digital initial Spectrul semnalului digital esantionat Spectrul semnalului decimat 53 Esantionarea spectrului unui semnal discret de durata finita Discretizarea semnalelor a fost impusa de utilizarea calculatoarelor numerice (care nu pot prelucra decat marimi discrete in timp/frecventa) Prelucrarea in domeniul spectrului este de multe ori mai simpla decat in domeniul timp Se pune intrebarea cu ce frecventa sa esantionam spectrul X(Ω) unui semnal discret x[n]                                  s s ssX X X k X k k kk      54 27 12      n       sN s sN 1 1 11              ;x n X X x n n                     N ss s 2  s N Fie semnalul x[n] cu suportul 0,M-1 Semnalul obtinut prin esantionarea spectrului semnalului x[n] este o extensie prin repetare cu perioada N a semnalului original x[n] N x n x n kN      2k    Conditia de reconstructie: x[n] sa fie de durata finita, MN 55 56 28 Se esantioneaza spectrul semnalului de durata finita x [n] (zero   pentru n M-1)  semnal periodic , perioada  nx~ 2 /sN Conditia de reconstructie : MN, nu sunt suprapuneri ale grupurilor temporale din semnalul rezultat •Reconstructia se face prin inmultirea semnalului cu fereastra temporala rectangulara: 2  , 0 1nN    x n x n x n w n w nN;             r r r  0, in rest  Daca spectrul se esantioneaza prea rar, rezulta N<M  apare suprapunerea grupurilor temporale; deci erori de aliere Semnalul x[n] nu mai poate fi reconstruit din spectrul esantionat 57 Masuri practice la esantionarea semnalelor analogice •Nu se cunoaste largimea benzii unui semnal de banda limitata ce urmeaza a fi esantionat •Semnalul poate avea componente spectrala de frecvente mari, neinteresante in aplicatia considerata Ele pot fi cauzate si de zgomot •Apare astfel riscul erorilor de aliere Acestea se evita folosind un filtru trece-jos numit si filtru anti-aliere prevazut in structura lantului de prelucrare a semnalului, inainte de circuitului de esantionare si memorare 58 29 •Sisteme de telefonie numerica (semnal audio) : •Frecventa maxima din spectru: f=3 4 kHz M •Frecventa Nyquist de esantionare: f=6 8 kHz s •Frecventa standard de esantionare: f=8 kHz s •Sisteme de televiziune (semnalul video): •Frecventa maxima din spectru: f=5 MHz M •Frecventa Nyquist de esantionare: f=10 MHz s •Frecventa standard de esantionare: f=18 MHz s 59 Semnal vocal fara aliere Semnal vocal cu aliere Semnal muzical fara aliere Semnal muzical cu aliere 60 30 Esantionarea semnalelor trece-banda Semnale de tip "trece jos" - spectrul concentrat in benzi care includ frecventa nula Semnale de tip "trece banda" - au suportul spectrului de forma supp X , ,                mmMM      Reconstructia perfecta a unui semnal trece banda esantionat ideal se poate realiza pe baza teoremei WKS: ω2ω sM Uneori semnalele trece banda pot fi reconstruite din esantioanele lor chiar daca s-a folosit o frecventa de esantionare mai mica decat frecventa Nyquist 61 Cazul semnalelor trece banda de banda ingusta   mM Semnal trece banda de banda ingusta : 1 m Suportul spectrului unui semnal trece banda de banda ingusta esantionat ideal este de forma: supp X , ,n n n n                      mmMMs s s s s     n Aici =  sM Reconstructie perfecta: filtrare trece-banda ideala chiar daca s-a folosit o frecventa de esantionare inferioara frecventei Nyquist 62 31 Filtrul trece-banda ideal 1,      12cc H p p             21ccTB 0, in rest  21sin sincctt   BPht tt 63 Reconstructie perfecta : , , ,k k l l k l                      mmMMs s s s     m, M] cu ls M,-m] cu ks Deplasare in intervalul [ Deplasare in intervalul [- Pentru l=0:                   , ,mmMMssk k k ,         22mk  MmMs      s 11kkk       MMs  Daca exista valori intregi ale lui k pentru care aceasta conditie este satisfacuta, atunci exista valori ale frecventei de esantionare inferioare frecventei Nyquist pentru care semnalele trece banda de banda ingusta pot fi reconstruite in urma esantionarii ideale 64 32 Solutia din multimea numerelor intregi a dublei inecuatii obtinute este m   00kn   mM frecventa de esantionare va apartine unor intervale de forma:   22mM   cu 1, ,kn0,  1kk  Exemplu m 10    04n 8msi atunci  Valorile M    mM admisibile pentru k sunt 1,2,3, 4 Intervalele corespunzatoare pentru frecventa de esantionare: 4 5 , 5 33 6 66 , 8 10 , 16 20 ,                        s  65 33 